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iy~j . Two, non-interacting two-level atoms immersed in a common bath can become 

mutually entangled when evolving with a Markovian, completely positive dynamics. 
For an environment made of external quantum fields, this phenomenon can be studied 
in detail: one finds that entanglement production can be controlled by varying the 
bath temperature and the distance between the atoms. Remarkably, in certain circum- 
stances, the quantum correlations can persist in the asymptotic long-time regime. 
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1 Introduction 



Independent atoms immersed in external quantum fields and weakly coupled to them can 
be viewed as open systems, i.e. as subsystems in interaction with an environment PQ-jl]. 
The atoms can be usually treated in a non-relativistic approximation, as independent n-level 
systems, with negligible size, while the environment is described by a set of quantum fields 
(e.g. the electromagnetic field) in a given quantum state, typically either a temperature 
state or simply the vacuum state. The interaction of the fields with the atoms is taken to 
be of dipole type, a well justified approximation within the weak coupling assumption j^J. 

Even in this simplified setting, that ignores all intricancies related to the internal atom 
structure and the full coupling with the electromagnetic field, the model is of great relevance 
both theoretically and phenomenologically 0-0 : indeed, with suitable adaptations, it is 
able to capture the main features of the dynamics of very different physical systems, like 
ions in traps, atoms in optical cavities and fibers, impurities in phonon fields. 

Despite this ample range of possible applications and the attention devoted to them in the 
recent literature, no particular care has often been taken in the derivation of an acceptable 
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sub dynamics for the atoms. As a result, time evolutions that are not even positive, have 
been adopted in order to describe their physical properties. 

On the other hand, a mathematically sound and physically consistent time evolution for 
the atom subsystem can be obtained using the weak coupling limit procedure the resulting 
subdynamics is described by a one parameter (= time) family of completely positive maps 
that form a quantum dynamical semigroup. 

In the following we shall outline such a derivation, and apply the resulting dynamics to the 
study of the evolution of a system composed by two, independent atoms. For simplicity we 
shall restrict the attention to two- level atoms in interaction with a collection of independent, 
free, massless scalar fields in 3 + 1 space-time dimensions, assumed to be in a state at 
temperature T — 1/(3. 

The interaction with an environment usually leads to decoherence and noise, typical mixing 
enhancing phenomena. Therefore, one generally expects that when a bipartite system is 
immersed in an environment, quantum correlations that might have been created before by 
a direct interaction between the two subsystems actually disappear. 

However, an external environment can also provide an indirect interaction between other- 
wise totally decoupled subsystems and therefore a mean to correlate them This 
phenomenon has first been established in exactly solvable models |0] : there, correlations be- 
tween the two subsystems take place during a short time transient phase, where the reduced 
dynamics of the subsystems contains memory effects. 

Remarkably, entanglement generation may also occur in the Markovian regime, through a 
purely noisy mechanism ITSt ITo] . It is precisely this situation that is relevant in the analysis 
of the dynamics of two independent atoms interacting with the same set of quantum fields. In 
the following, we shall study in detail the conditions that allow the two otherwise indepedent 
atoms to become initially entangled through the action of the environment, paying special 
attention to the external controllable parameters, the bath temperature and the spatial 
distance i between the two atoms. We shall see that for fixed, finite £, there is always a 
temperature below which entanglement generation occurs as soon as time starts to become 
nonzero. Remarkably, it is found that for vanishing i the entanglement thus generated 
persists even in the long-time asymptotic equilibrium state. 

2 Two Atom Master Equation 

We shall deal with a system composed by two, identical two-level atoms, that start inter- 
acting at time t = with a collection of independent, massless, scalar quantum fields at 
temperature T. We are not interested in the details of the atoms internal dynamics. We 
shall therefore model them, in a nonrelativistic way, as simple two-level systems, which can 
be fully described in terms of a two-dimensional Hilbert space. 

In absence of any interaction with the external fields, the single atom internal dynamics 
will be driven by a 2 x 2 hamiltonian matrix, that in a given basis can be taken to assume 
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the most general form: ■ 3 = ^J2i=i n i cr i , where <7j, 2 = 1,2,3 are the Pauli matrices, 
rij, i — 1,2,3 are the components of a unit vector, while u; represents the gap between the 
two energy eigenvalues. Then, the atom Hamiltonian Hg is the sum of two such terms: 

Hs = Hjp + HV, Hf = £2>a<°> , a = 1,2, (1) 

where = a, ® 1 and af^ = 1 <8> cr, are the basis operators pertaining to the two different 
atoms. 

As mentioned in the Introduction, the interaction of the atoms with the external fields is 
assumed to be weak; it can then be described by an hamiltonian H' that is linear in both 
atom and field variables: 

ff , = E(^ 1) ®$<[/ (1) ]+^ 2) ®$i[/ (2) ]) • (2) 
i=i 

The operators &i(t,x) represent the external quantum fields, taken to be spinless and mass- 
less for simplicity. They evolve in time as free relativistic fields with a standard Hamiltonian 
H$. The atoms are assumed to have a spatial extension described by the two functions 
f^ a \x), a = 1,2, taken to have a common profile f(x). To be more specific, we shall choose 
for the atoms a spherically symmetric shape of infinitesimal size e: 



(f/2) 
7T 2 [|£|2 + (e/2) 2 }' 



m = - ,, . • (3) 



Further, without loss of generality, the first atom can be positioned at the origin of the 
reference frame, so that one can assume f^\x) = f(x), while the second is displaced by 
an amount £ with respect to it, and therefore: f^(x) = f(x + £). Since the atom-field 
interaction takes place on the whole region occupied by the atoms, the field operators entering 
the interaction Hamiltonian above are smeared over the atom size: 

^[f {a) ] = J d 3 xf^(x)U0,x) , a= 1,2 . (4) 

The total Hamiltonian H describing the complete system, the two atoms together with 
the external fields $j, can thus be written as 

H = H s + + XH' = H + XH' , (5) 

with A a small coupling constant. It generates the evolution in t of the corresponding total 
density matrix p tot , dtPtot(t) — —i[H, ptot(^)], starting at t — from the initial configuration: 
Ptot(O); we shall assume the atom and the fields to be initially prepared in an uncorrelated 
state, with the fields in the temperature state pp and the atoms in a generic initial state 
p(0), so that ptot(O) = p(0) <g> pp. 



3 



At this point, being interested in studying the dynamics of the two atoms, one conveniently 
integrates over the unobserved field degrees of freedom and concentrate on the analysis of 
the reduced time evolution, formally given by the transformation map: p(0) ^ p(t) = 
Tr$[p t ot(*)]- 

The derivation of a physically consistent master equation for the reduced density matrix 
p{t) = Tr$[p tot (t)] is notoriously tricky, requiring an a priori unambigous separation between 
subsytem and environment, besides a sufficiently weak interaction between the two 
Generally speaking, this distinction can be achieved when the correlations in the environment 
decay much faster than the characteristic evolution time of the subsystem alone, given by 
the inverse of its typical energy scale. In such a case, in the limit of weak couplings, the 
changes in the evolution of the subsystem occur on time scales that are very long, so large 
that the details of the internal environment dynamics result irrelevant. This is precisely the 
situation that occurs for the system under study: indeed, in typical instances, the differences 
between the atomic internal energy levels result much smaller than the field correlation decay 
constants so that a clear distinction between subsystem and environment is authomatically 
achieved. 

In practice, the dynamics of the reduced system is obtained by suitably rescaling the time 
variable, t — > t/X 2 and then taking the limit A — * 0, following the mathematically precise 
procedure of the weak coupling limit |H],[I]-jSl- The reduced density matrix p(t) is then found 
to obey the following evolution equation: 

^ = L Hs [p(t)] + > L Hs [p] -i[H s , p] , (6) 

where 



1 f T 

V*[- \ = - lim 

T — >CO 

and 



Z>«[-] = - lim - / dsU{-s)VU{s)[-\ , U(s) = e sL »s, (7) 

T-*oo 1 Jo 



V[p\ = J™ dt Tr([e iHot H'e- lHot } [H',p® Pf) ]\) . (8) 

For the case at hand, the integrals in (J7J) and (JBJ) can be explicitly computed and the master 
equation for p{t) written down without any ambiguity. It takes a Kossakowski-Lindblad 
form [13 HE! 

cW) --i[H eS ,p{t)]+C\p{t)\ , (9) 



dt 



with 



and 



2 3 



H« = Hs-\Y. E H\f ] 4 a) °f , (10) 

a,P=li,j=l 



1 



C[P] = ;EE Clf^afpaf-aPafp-paPaf] . (11) 

A a,P=l i,j=l 
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The coefficients of the Kossakowski matrix C\"^ and of the effective Hamiltonian H e g are 
determined by the field correlation functions in the thermal state pp. 

G \f\t -t') = J d 3 xd 3 yf^(x) f^(y) x)^(t', y)) , (12) 

through their Fourier, 

G\f\z)= r dte^G\f\t) , (13) 



and Hilbert transform, 

)C\f\z) = r dtsign(t)e izt G\f\t) = - f°° dw , (14) 

3 J -oo 3 Til J-00 W - Z 

respectively (P indicates principle value). More specifically, one finds: 



c\f ] = E E^(^)4?W- 5) , (is) 

£=+,-,0 k,l=l 



and similarly for H^j, with Q^\^uj) replaced by /Cjy (£a;), where 



4>-f = rii rij , iplf = -(Sy - rii rij ± it ijk n k ) . (16) 



are the components of auxiliary three-dimensional tensors. 1 Being the sum of three positive 
terms, the matrix C-?® turns out to be positive, so that the dynamical semigroup generated 
by is composed by completely positive maps: this is the result of adopting a mathemat- 
ically well-defined formalism jHj. On the other hand, let us remark that direct use of the 
standard second order perturbative approximation {e.g. see 0111) often leads to physically 
inconsistent results, giving a finite time evolution for p(t) that in general does not preserve 
the positivity of probabilities. 

The expressions in ()13|) and can be explicitly computed by noting that the fields are 
taken to be independent and assumed to obey a free evolution, so that: 

= ^IM^Mpp] = % G ( x - v) , ( 17 ) 

where G(x — y) is the standard four- dimensional Wightmann function for a single relativistic 
scalar field in a state at inverse temperature f3, formally given by: 



d^k 



G(x) = J 6(k°) 5(k 2 ) [(1 + M(k )) e- tk ' x + Af(k°) e %k ' r 

where 



1 We omit the details of the derivation and refer to the Appendix of JHj for an outline of the needed 
techniques. 
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Note that, as a result, the correlations in (fT2*j) involve f(k) = J d 3 x e lk ' x f(x), namely the 
Fourier transform of the shape function f(x) in 0; it can be easily computed to be f(k ) = 
e -\k\e/2_ B e j n g a function of the modulus \k\ only, this contribution can be conveniently 
attached to the definition of the Wightmann function G(x), so that the integrand in (fTBj) 
gets an extra e~ £k ° overall factor. This damping term assures now the convergence of the 
integral ()18|) and corresponds to the usual is prescription for the Wightmann function; in 
the present setting, it 0X1 S6 clS db remnant of the (infinitesimal) size of the atoms. 

The behaviour of the field correlations in (|T2*|) is also crucial for assuring the convergence 
of the evolution equation of the reduced density matrix p(t) to the limit ©; indeed, one 
shows jHj that such limit exists only when the combination |G^°^(t)|(l + t) v is integrable on 
the positive half real line, for some rj > 0. In the case of massless fields considered here, this 
condition is assured by the 1/t 2 fall off at infinity of the Wightmann function. 

Using (fTTj) and (fTHj) . the Fourier transform in (fT3~j) can now be explicitly evaluated (in the 
limit of vanishing e) : 



with: 



gM( z ) = g<V)( z ) 



(20) 



z 



2tt 1 - e-P* ' 
1 z sm(£z) 



2tt 1 - e-P* iz 



(21) 



where i denotes the modulus of the displacement vector £; then, recalling (pPf]) . for the Hilbert 
transform one similarly finds: 



J m J-oo w — z 

With these results, The Kossakowski matrix can finally be written as: 



(22) 



r>( 12 ) 



Clj 2) = A 5^ - iB e ijk rik + Crii n j , 
C\f ] = A' - iB' e ijk n k + C rij 



3 ' 



(23) 



where the quantities A, B, C, A', B' and C depend on the system frequency u, the inverse 
temperature (3 and the separation i between the two atoms: 



.4 
B 
C 



Air 
Air 

UJ 

Air 



1 + e-P" 



■ 2 

.Qui 



-f3u 



1 + e 



-0LJ- 



1 - e~^. 



A' 



B' 



C 



Air 



1 + e _/Jw 



-put 



sin (ojl 



Air lo£ 



UJ 

Air 



1 + e 



sin wt 



_(3uj 1 - e~P" uj£ 



(24) 
(25) 
(26) 
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On the other hand, the effective Hamiltonian H c r naturally splits into two contributions: 

~ (12) 

H e ff = Hs + i/eff • The fi rs * one ^ as the same form of the starting system Hamiltonian in 
(JTJ), but with a redefined frequency: 

Co = 00 + i[K, {ll) (-uo) - /C (11) (oo)} , (27) 

while the second one corresponds to an environment generated direct coupling among the 
two atoms: 



H^ = - l -j: ([IC^(-oo) + IC^(oo)]5., 

z i,i=l 



IJ 



+ [2/C (12) (0) - /C (12) (-^) - /C (12) M] mrij) (n ® a, . (28) 

These results for the Hamiltonian contributions require some further comments. Recalling 
(JUJ), the definition of /C^ 11 ^^) in (J2~2~j) can be split as (similar results hold also for K,^ (z)): 



^ u \z) -- ' 



2txH 



oo 



P / ^ (29) 
w — 2 



+P / diw 



..... , , ■ (30) 

1 — e^ w \ w + z w — z ' 

into a vacuum and a temperature-dependent piece. Although not expressible in terms of 
elementary functions, the temperature dependent second term is a finite, odd function of 
z, vanishing as (3 becomes large, i.e. in the limit of zero temperature (and as such, it does 
not contribute to H^ s in (|28j)). The first contribution in (j2*D*j) is however divergent. As a 
consequence, despite some cancellations that occur in (j2*7|) and (|28j) . the effective Hamiltonian 
H e ff turns out in general to be infinite, and its definition requires the introduction of a suitable 
cutoff and a renormalization procedure. 

This is a well known fact, and has nothing to do with the weak-coupling assumptions used 
in deriving the master equation. Rather, the appearance of the divergences is due to the 
non-relativistic treatment of the two-level atoms, while any sensible calculation of energy 
shifts would have required the use of quantum field theory techniques 

In our quantum mechanical setting, the procedure needed to make H e R well defined is 
therefore clear: perform a suitable temperature independent subtraction, so that the ex- 
pressions in (|27|) and ()28)1 reproduce the correct quantum field theory result, obtained by 
considering the fields in the vacuum state. 

In the following, we shall be interested in analyzing the temperature dependent effects 
described by the master equation 0; all standard, vacuum generated Hamiltonian contri- 
butions will be therefore ignored. 



3 Environment Induced Entanglement Generation 

Using the explicit form for the master equation derived in the previous Section, one can 
now investigate whether the thermal bath made of free fields can actually entangle the two 
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indepedent atoms. Since we are dealing with a couple of two-level systems, this can be 
achieved with the help of the partial transposition criterion fOJ|2n|: a two-atom state p(t) 
results entangled at time t if and only if the operation of partial transposition does not 
preserve its positivity. 

We shall first consider the possibility of entanglement creation at the beginning of the 
evolution. Without loss of generality one can limit the considerations to pure, separable 
initial states, and therefore take: 

p(0) = \<p)(<p\®\1>)(iJ>\ ; (31) 

indeed, if the environment is unable to create entanglement out of pure states, it will certainly 
not entangle their mixtures. Then, let us examine the behavior in a neighborhood of t = 
of the quantity 

Q(t) = (Xl m \X) , (32) 
where the tilde signifies partial transposition, e.g. with respect to the second factor, and \x) 
is any 4-dimensional vector. The two atoms, initially prepared in a state p(0) = p(0) as in 
(J2J), will surely become entangled if there exists a suitable vector \x), such that: i) <2(0) = 
and ii) <9t<2(0) < 0. In fact, when d t Q(0) > for all choices of the initial state p(0) and probe 
vector \x), entanglement can not be generated by the environment, since p remains positive. 
Clearly, the vector \x) need be chosen entangled, since otherwise Q is never negative. 

Note that study of the behavior of the quantity Q(t) near t = allows a very explicit 
analysis of entanglement generation. Indeed, d t Q(0) can be easily computed through the 
time derivative d t p(0), that in turn can be obtained by taking the partial transposition of 
the r.h.s. of Q (with if e g set to zero, as explained above). In this way, one can construct a 
test of entanglement creation, valid for any probe vector \x). 

In order to show this, consider first the orthonormal basis {\<p), \<f)},{\ip), obtained 
by augmenting with the two states \(p) and the ones that define p(0) in 1)31)). They can 
be both unitarily rotated to the standard basis {|— ), |+)} of 03: 

w) = u\~) m = u\+) , 

W) = V\-) $) = V\+). (33) 

Similarly, the unitary transformations U and V induce orthogonal transformations U and V, 
respectively, on the Pauli matrices: 

3 3 
rfaiU = £ UiM , V^tnV = £ Vijffj . (34) 
3=1 i=i 

Direct computation then shows that d t Q(0) can be written as a quadratic form in the in- 
dependent components of the probe vector \x), with coefficients that involve the four 3x3 
matrices C^ ll \ C^ 22 \ C^ l2 \ that form the Kossakowski coefficients given in (jl5j) . As a 
consequence, vectors \x) exist making this form negative, i.e. d t Q(0) < 0, if and only if its 
corresponding discriminant is negative; explicitly: 

{u\C {11) \u) (v\(C {22) ) T \v) < \(u\1le(C {l2) )\v)\ 2 , (35) 
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where T means matrix transposition. The three-dimensional vectors \u) and \v) contain 
the information about the starting factorized state (|31|) : their components can be in fact 
expressed as: 

3 3 

(+\aj\-) , Vi = £ V« <-N+> • (36) 



U: 



Therefore, the external quantum fields will be able to entangle the two atoms evolving with 
the Markovian dynamics generated by Q and characterized by the Kossakowski matrix (jl5j) . 
if there exists an initial state \<p){<p\ <8> \ip)(ip\ ) or equivalently orthogonal transformations U 
and V, for which the inequality (J3*5)l is satisfied. 

That this is indeed the case for the matrices in (J23)l can be easily shown. First note 

that, without loss of generality, the unit vector n that defines the internal atom Hamiltonian 
in (|TJ) can be oriented along the third axis. Consider then the initial state p(0) = |— )(— | <8> 
| +)(+|, constructed out of the eigenstates of the single atom Hamiltonian. Recalling the 
definitions (|3Tjl and (|33|). one then finds U = 1, so that the three-dimensional vector \u) has 
components U{ = {1, —i, 0}, and further Vi = Ui. Then, the inequality reduces to: 

rf + S »>l, fl 3 * = i^±, s=2^. (37) 

A 1 + e-P" cu£ K J 

Although both R and S take values in the interval [0, 1], one can easily make the sum of their 
squares to exceed unity by adjusting the inverse temperature f3 and the atom separation £. 
In particular, for a given, finite separation £, one can always find a temperature below which 
the inequality in (|3"7jl is satisfied and therefore entanglement created. The case of vanishing 
separation is even more striking, since the inequality (J3"7)l reduce to R > 0, which is always 
satisfied, except in the limit of an infinite bath temperature. 

Analogous results hold also in the case of a zero temperature bath. Indeed, in this case 
R = 1, and entanglement is generated for any finite separation of the two atoms (similar 
conclusions have also been reported before in Ref. \21\). 

The role played by the two considered control parameters, the bath temperature and the 
atom separation, in triggering entanglement creation is now apparent. The temperature of 
the external environment determines the amount of noise that is induced in the dynamics of 
the two independent atoms. Noise is known to reduce quantum correlations, and indeed, the 
higher the temperature, the less effective is the entanglement power of the bath. Environment 
induced entanglement generation is nevertheless a robust phenomenon: it always occurs 
except in the limit of an infinitely large temperature. 

A similar role is played by the second control parameter, the spacial atom separation: 
entanglement enhancement is more effective the closer the two atoms are, and turns out to 
be impossible only for an infinitely large separation. The interplay between the effects of 
these two control parameters is neatly summarized by the inequality (j3"7l) . 

To our knowledge, this is the first instance of open quantum system control through the 
bath parameters, and not via the subsystem Hamiltonian fl^\ 123"] . This approach, which 
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might prove very fruitful in the field of quantum information, is still in a very preliminary 
stage; further developments are presently under study and will be reported elsewhere. 

Finally, note that the above choice for p(0) is not restrictive: one can always use the trans- 
formations in (J3~3~j) to map it to the generic initial state (JUJ), at the expences of using suitably 
rotated \u) and \v), as given in (jHEJ)- As a result, the expression of the entanglement test in 
(jHSJ) becomes more complicated than in (jHZJ), but the final conclusions remain unchanged. 



4 Asymptotic Entanglement 

On the basis of the anlysis presented in the previous Section two atoms, initially prepared in a 
separate state, will generically become entangled as a result of their independent interaction 
with a bath made of thermal quantum fields: quantum correlations are generated among the 
two atoms as soon as t > 0. The test in (J2HJ), on which this conclusion is based, is however 
unable to determine the fate of this quantum correlations, as time becomes large. 

On general grounds, one expects that the effects of decoherence and dissipation that 
counteract entanglement production be dominant at large times, so that no entanglement 
is left in the end. This is precisely the conclusion that one obtains by a careful analysis of 
the structure of the dynamics generated by the master equation in ©, with Kossakowski 
coefficients as in (}23|) - (j26|) : the asymptotic equilibrium state for the two atoms turns out to 
be always separate for £ finite (a detailed account of this result is beyond the scope of the 
present work, and will be reported elsewhere). 

However, the case of a vanishing atom separation is again special and deserves a closer 
examination. Note that in such a situation, the matrices in (|23|) become all equal C^- 11 ' = 
q{22) _ (j{i2) _ £7(21) _ Thig particular choice for the Kossakowski matrix is also adopted in 
the description of the phenomenon of resonance fluorescence |24[ I7]. Therefore, the discussion 
of the vanishing i limit is of relevance also from the phenomenological point of view. 

The presence of an equilibrium state p°° can be in general determined by setting to zero the 
r.h.s. of the evolution equation Q. As previously explained, we shall ignore the Hamiltonian 
piece since it can not give rise to temperature dependent entanglement phenomena, and 
concentrate on the study of the effects induced by the dissipative part; the equilibrium 
condition reduces then to C[p°°] = 0. Direct computation leads to the following result: 

1 r 3 3 i 

5 1 — a n iiX ® <7i + °i ® 1) + X! {bdij + CTliTlj) <Ji ® Oj 



=1 



(38) 



with 

a= jftf( T + 3 ^ 6= ^H^ r + 3 )' c = Ra - ( 39 ) 

Here, R = B/A is the temperature dependent ratio already introduced in (|3*7j). while the 
parameter r = \ Z)E=i Tr[p(0) (crjCgxjj)] contains the dependence on the initial state (positivity 
of p(0) requires —3 < r < 1): the equilibrium state is therefore not unique. In spite of this, 
remarkably, the asymptotic state (jSH|) turns out to be still entangled. 
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To explicitly show this, one can as before act with the operation of partial transposition 
on p°° to see whether negative eigenvalues are present. Alternatively, one can resort to 
one of the available entanglement measures and concurrence appears here to be the more 
appropriate: its value C[p] ranges from zero, for separable states, to one, for fully entangled 
states, like the Bell states (2i]-j2Zl- In the case of the state p°° above, one finds 



C[p°°] = max< 



(3 - R 2 ) 
2(3 + B?) 



5R 2 -3 
3-R 2 ~ T . 



. (40) 



This expression is indeed nonvanishing, provided we start with an initial state p(0) for which 

T< Y^W (41) 

The concurrence is therefore a linearly decreasing function of r, starting from its maximum 
C[p°°} = 1 for r = -3 and reaching zero at r = (5i? 2 - 3)/(3 - R 2 ). 

This result is remarkable, since it implies that the dynamics in (JSJ) not only can initially 
generate entanglement: it can continue to enhance it even in the asymptotic long time 
regime. 2 In other terms, put the two atoms in the same place and prepare them at t — in 
a separable state; then, provided the condition (j41j) is satisfied, their long time equilibrium 
state will turn out to be entangled. 

The simplest example of a separable state is provided by the direct product of pure states 
as in (|31|). When \ip) and are orthogonal, so that r = —1, one explicitly finds: C[p°°] = 
max{2i? 2 /(3 + R 2 ), 0}. Notice that C[p°°] reaches its maximum value of 1/2 when R = 1, 
i.e. at zero temperature, while it vanishes when the temperature becomes infinitely large, 
i.e. R = 0. As already remarked, this has to be expected, since in this case the decoherence 
effects of the bath become dominant. 



References 

[1] V. Gorini, A. Frigerio, M. Verri, A. Kossakowski and E.G.C. Sudarshan, Rep. Math. 
Rhys. 13 (1978) 149 

[2] H. Spohn, Rev. Mod. Phys. 52 (1980), 569 

[3] R. Alicki and K. Lendi, Quantum Dynamical Semigroups and Applications, Lect. Notes 
Phys. 286, ( Springer- Verlag, Berlin, 1987) 

[4] H.-P. Breuer and F. Petruccione, The Theory of Open Quantum Systems (Oxford Uni- 
versity Press, Oxford, 2002) 

2 Similar results hold in the case of uniformily accelerating atoms ^Bj; in that case, entanglement gener- 
ation is controlled by the value of the proper acceleration. 



11 



[5] P.W. Milonni, The Quantum Vacuum: An Introduction to Quantum Electrodynamics, 
(Academic Press, San Diego, 1994) 

[6] Z. Ficek and R. Tanas, Phys. Rep. 372 (2002) 369 

[7] R.R. Puri, Mathematical Methods of Quantum Optics, (Springer, Berlin, 2001) 

[8] E.B. Davies, Comm. Math. Phys. 39 (1974) 91; Math. Ann. 219 (1976) 147 

[9] D. Braun, Phys. Rev. Lett. 89 (2002) 277901 

[10] M.S. Kim et al, Phys. Rev. A 65 (2002) 040101(R) 

[11] S. Schneider and G.J. Milburn, Phys. Rev. A 65 (2002) 042107 

[12] A.M. Basharov, J. Exp. Theor. Phys. 94 (2002) 1070 

[13] L. Jakobczyk, J. Phys. A 35 (2002) 6383 

[14] B. Reznik, Found. Phys. 33 (2003) 167 

[15] F. Benatti, R. Floreanini and M. Piani, Phys. Rev. Lett. 91 (2003) 070402 

[16] F. Benatti and R. Floreanini, Phys. Rev. A 70 (2004) 012112 

[17] V. Gorini, A. Kossakowski and E.C.G. Sudarshan, J. Math. Phys. 17 (1976), 821; 

[18] G. Lindblad, Comm. Math. Phys. 48 (1976) 119 

[19] A. Peres, Phys. Rev. Lett. 77 (1996) 1413 

[20] M. Horodecki, P. Horodecki and R. Horodecki, Phys. Lett. A 223 (1996) 1 

[21] Z. Ficek and R. Tanas, J. Mod. Opt. 50 (2003) 2765 

[22] S. Lloyd and L. Viola, Phys. Rev. A 65 (2002) 010101 

[23] C. Altafini, J. Math. Phys. 44 (2003) 2357 

[24] G.S. Agarwal, A.C. Brown, L.M. Narducci and G. Vetri, Phys. Rev. A 15 (1977) 1613 

[25] S. Hill and W.K. Wootters, Phys. Rev. Lett. 78 (1997) 5022 

[26] W.K. Wootters, Phys. Rev. Lett. 80 (1998) 2245 

[27] W.K. Wootters, Quantum Inf. Comp. 1 (2001) 27 



12 



